I- 



(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property Organization 

International Bureau 

(43) International Publication Date 
30 October 2003 (30.10.2003) 




(10) International Publication Number 

PCT WO 03/089954 Al 



(51) International Patent Classification 7 : G01S5/14, 
H04Q7/20 

(21) International Application Number: PCT/US03/11475 

(22) International Filing Date: 1 1 April 2003 (1 1 .04.2003) 

(25) Filing Language: English 

(26) Publication Language: English 



(30) Priority Data: 
10/126,035 



18 April 2002 (18.04.2002) US 



(71) Applicant: QUALCOMM, INCORPORATED 

[US/US]; 5775 Morehouse Drive, San Diego, CA 92121 
(US). 

(72) Inventor: GAAL, Peter; 95 15 Easter Way, #6, San Diego, 
CA 92121 (US). 

(74) Agents: WADSWORTH, Philip R. et al.; 5775 More- 
house Drive, San Diego, CA 92121 (US). 



(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, FI, GB, GD, GE, GH, 
GM, IIR, HU, ID, IL, IN, IS, JP, KB, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, 
MX, MZ, NI, NO, NZ, OM, PH, PL, PT, RO, RU, SC, SD, 
SB, SG, SK, SL, TJ, TM, TN, TR, TT, TZ, UA, UG, UZ, 
VC, VN, YU, ZA, ZM, ZW. 

(84) Designated States (regional): ARTPO patent (GH, GM, 
KE, LS, MW, MZ SD, SL, SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE, BG, CH, CY, CZ, DE, DK, EE, 
ES, FI, FR, GB, GR, HU, IE, IT, LU, MC, NL, PT, RO, 
SB, SI, SK, TR), OAPI patent (BF, BJ, CF, CG, CI, CM, 
GA, GN, GQ, GW, ML MR, NE, SN, TD, TG). 

Published: 

— with international search report 

— before the expiration of the time limit for amending the 
claims and to be republished in the event of receipt of 
amendments 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations'' appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



(54) Title: INTEGRITY MONITORING IN A POSITION LOCATION SYSTEM UTILIZING KNOWLEDGE OF LOCAL TO- 
POGRAPHY 



00 

rn 
O 



DERIVING PLURALITY 
OF POSITION 
ESTIMATES, EACH 
FROM A DIFFERENT 
MEASUREMENT 
SUBSET 



SELECTING POSITION 
ESTIMATE BASED ON 
LOCAL TOPOGRAPHY 



1100 



.1102 



(57) Abstract: A system for and method of deriving a position estimate of a 
receiver from a plurality of parameter measurements is provided. Each mea- 
surement is derived from a correlation function, which in rum is derived by cor- 
relating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver. A plurality of position esti- 
mates are derived, each from a different subset of the plurality of measure- 
ments. Based on local topography, one of the position estimates is selected as 
the best estimate of the position of the receiver. 
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INTEGRITY MONITORING IN A POSITION LOCATION SYSTEM 
UTILIZING KNOWLEDGE OF LOCAL TOPOGRAPHY 



1. Field of the Invention 

[0001] This invention relates to the field of integrity monitoring of measurements taken or 
position estimates derived in a position location system, and, more specifically, to performing 
integrity monitoring utilizing knowledge of local topography to account for factors such as lack 
of line-of-sight and multi-path. 

2. Related Art 

[0002] The Global Positioning System (GPS) is a collection of satellites each of which 
travels in a precise orbit above the earth's surface. Each satellite transmits a signal modulated 
with a pseudo-noise (PN) code unique to the satellite. Each PN code comprises a predetermined 
number of chips. A GPS receiver receives a composite signal comprising a mixture of signals 
from each of the satellites that are visible to the receiver. A signal detector in the receiver 
detects a transmission from a particular satellite by determining the degree of correlation 
between the received signal and shifted versions of the PN code for that satellite. If a peak of 
sufficient quality in the correlation value for one of the shift offsets is detected, the receiver is 
considered to have detected the transmission from the satellite. 

[0003] The receiver estimates its location by detecting transmissions from at least four of 
the satellites. For each detected transmission, the receiver uses the shift in the PN code to 
estimate the delay (in terms of chips or fractions of chips) between time of transmission and 
time-of-arrival. Given the known propagation speed of the transmitted signal, the receiver 
estimates the distance between itself and the satellite. This estimated distance defines a sphere 
around the satellite. The receiver knows the precise orbits and positions of each of the satellites, 
and continuously receives updates to these orbits and positions. From this information, the 
receiver is able to determine its position (and the current time) from the point where the spheres 
for the four satellites intersect. 

[0004] The FCC has mandated that subscriber stations, including but not limited to 
mobile stations, in wireless communications systems be capable of estimating their locations in 
order to promote rapid responses to 91 1 and other emergency calls. In response to this mandate, 
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efforts are underway to equip subscriber stations with the means to estimate their locations from 
GPS satellite transmissions. Moreover, since base stations or sectors in wireless 
communications systems transmit pilot signals modulated with unique PN codes, these efforts 
also include allowing subscriber stations to estimate their locations from the transmissions of 
multiple base stations or sectors, or combinations of base stations or sectors and GPS satellites. 
The use of base stations or sectors to estimate position is referred to as Advanced Forward Link 
Trilateration (AFLT); and the use of both GPS satellites and base stations or sectors at the same 
time is referred to as hybrid position location. 

[0005] Referring to Figure 1, a signal detector within subscriber station 100 attempts to 
detect the transmission of a reference source, whether it be a satellite, base station or sector, 
from the peak of a correlation function which is derived by multiplying the received signal 
(which is typically a composite signal comprising a mixture of the transmissions from multiple 
reference sources, such as base stations/sectors 102a, 102b, 102c, and/or satellites 106a, 106b, 
106c) with shifted versions of the PN code for the reference source within a range defined by a 
predetermined search window, and then, for each shifted PN code, adding the multiplied values 
over a predetermined integration time to achieve a value representative of the degree of 
correlation between the received signal and the shifted PN code. If a peak is detected, the signal 
detector may then derive measurements of one or more parameters, such as time-of-arrival, 
from the peak. Once four or more such measurements have been obtained, an estimate of the 
position of the subscriber station 100 may be obtained 

[0006] At least four measurements are required to determine a position estimate. If more 
than four measurements are obtained, the position estimate is considered to be over-determined. 
In an over-determined situation, it is possible to derive multiple position estimates, each from 
different subsets of the measurements. Further, if a measurement subset itself is over- 
determined, then it is possible to check that subset for self-consistency. In a process known as 
integrity monitoring, a quantitative measurement of the consistency between measurements in a 
subset is determined The quantitative measurement may be associated with the position 
estimate derived from that subset and be referred to as the 'cost' of the position estimate. These 
costs may be useful in selecting one of the estimates as the best estimate. 
[0007] This process may be further explained with reference to Figure 2. There, it is 
assumed that measurements have been obtained from reference sources 202a, 202b, 202c, 202d, 
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202e, 202f and that three position estimates, pei, pez, and pe 3 , have been obtained from different 
subsets of the measurements. In particular, as illustrated, pei is assumed to be derived from 
measurements taken from reference sources 202a, 202b, 202c, 202d and 202e; pe 2 is assumed to 
be derived from measurements taken from reference sources 202a, 202b, 202c, 202d, and 202f; 
and pe3 is assumed to be derived from measurements taken from reference sources 202b, 202c, 
202d, 202e, and 202f. 

[0008] In one example, the cost is the RMS error residual of the measurements used in 
deriving the position estimate. In this example, the cost associated with the first estimate, pei, 
may be expressed as follows: 

cost, =p„ -d aptl ) 2 +{r b -d bpcl ) 2 +(r c -d tpel ) 2 + (r rf -d d ^) 2 + (r„ -d e ^) 2 , 
(1) 

where r a , r b , r r , r d , and r e are the ranges measured between the subscriber station and 
reference sources 202a, 202b, 202c, 202d, 202e, respectively; and d a pet , d b ptX , d c pe) , d d pel , 
and d t pei are the computed distances between position estimate pei and reference sources 202a, 
202b, 202c, 202d, and 202e, respectively. 

[0009] Similarly, in this example, the cost associated with the second estimate, pe2, may 
be expressed as follows: 

cost 2 = J(r tt -d apt2 ) 2 + {r b -d bpt2 ) 2 +(r c -d cpe2 ) 2 + (r d -d dpt2 ) 2 +(r, -d /pt2 ) 2 , 
(2) 

where r f is the range measured between the subscriber station and reference source 202f; and 
d o pt2 , d bpe2 , d r pe2 , d d ^ , and rf /-pt2 are ^ computed distances between position estimate 
pe2 and reference sources 202a, 202b, 202c, 202d, and 202f, respectively. 
[0010] The cost associated with the third estimate, pe3, may be expressed as follows: 

cost, =J{r b -d bpe3 ) 2 +(r e -d cpc ,f + (r„ -d dpt3 ) 2 + (r, -d.^f + (r, -d,^) 2 , 
(3) 
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where d bpti , d cpe3 , d dptZ , d tpt) , and d f jxi are the computed distances between position 

estimate pe 3 and reference sources 202b, 202c, 202d, 202e, and 202f, respectively. 
[0011] In a conventional position location system, the position estimate corresponding to 
the lowest cost is selected as the best position estimate. However, there are many cases where 
the lowest cost estimate is not in fact the best estimate. Consider, for example, a scenario where 
more than one of the measurements used to compute the lowest cost estimate have been 
corrupted due to factors such as multi-path and lack of line of sight, but the measurements used 
to determine a higher cost estimate are not so corrupted. In this scenario, the higher cost 
estimate might in fact be the best position estimate. 

[0012] If the corrupted measurements could be identified, accuracy could be improved by 
excluding from consideration position estimates based on corrupted measurements. However, 
with the conventional approach, if less than a total of six measurements are available, or if more 
than one of the measurements have been corrupted, then it is generally impossible to identify 
which of the measurements have been corrupted. 

SUMMARY 

[0013] The invention provides a method of deriving a position estimate of a receiver from 
a plurality of parameter measurements. Each of the parameter measurements is derived from a 
correlation function, which in turn is derived by correlating a signal received at the receiver with 
an identification code uniquely identifying a transmitter visible to the receiver. 
[0014] A first embodiment comprises deriving a plurality of position estimates, each from 
a different subset of the plurality of measurements. One of the position estimates is then 
selected based on local topography. 

[0015] A second embodiment of the method comprises deriving a plurality of position 
estimates, each from a different subset of the plurality of measurements. The cost of each of the 
position estimates is then derived. The cost represents the degree of consistency of the 
measurements used in deriving the position estimate. Once obtained, the costs are then refined 
based on local topography. One of the position estimates is then selected based on the refined 
costs. 

[0016] A third embodiment of the method comprises using local topography to exclude 
any of the plurality of measurements which are unreliable. A plurality of position estimates are 
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then derived, each from a different subset of the remaining measurements. Costs for each of the 
position estimates are then determined, and one of the position estimates is then selected 
responsive to the costs. 

[0017] A fourth embodiment of the invention comprises deriving a plurality of position 
estimates, each from a different subset of the plurality of measurements. A reliability metric is 
then derived for each of the position estimates based on local topography. One of the position 
estimates is then selected based on the reliability metrics. 

[0018] In one implementation, the identification codes are pseudo noise (PN) codes, and 
the measurements are time-of-arrival measurements. In this implementation, a plurality of 
position estimates are derived, each from a different subset of the time-of-arrival measurements. 
[0019] Based on local topography, a reliability metric for each position estimate is derived 
using the following procedure. For each time-of-arrival measurement from which the position 
estimate is derived, based on local topography, the expected peak amplitude of the correlation 
function corresponding to a line of sight transmission is estimated. This value is then compared 
with the amplitude of the earliest non-sidelobe peak of the actual correlation function. The 
reliability metric is then adjusted based on the comparison. 

[0020] For example, if the actual peak amplitude exceeds the expected peak amplitude, a 
multi-path condition could be present, and the reliability metric decreased to reflect that. As 
another example, if the actual peak amplitude is less than expected, there might be unexpected 
obstacles in the path of the transmission from the transmitter, and the reliability metric 
decreased to reflect that This step is repeated for each of the measurements from which the 
position estimate is derived. The result is the reliability metric for the position estimate. 
[0021] Each of these methods may be tangibly embodied as a series of processor 
executable instructions stored a memory accessible by the processor. Furthermore, systems are 
possible comprising a processor and the foregoing memory, wherein the memory is accessible 
by the processor, and the processor is configured to access and execute the instructions stored 
thereon. 

[0022] Such a system may be used in a variety of environments or applications. For 
example, in one application, such a system forms a part of a subscriber station and is configured 
to locate the position of the subscriber station from one or more transmissions from base 
stations, sectors, GSP satellites, or any combination thereof. In another application, such a 
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system forms part of a position determination entity (PDE), which is configured to locate the 
position of a subscriber station from measurements obtained by the subscriber station. More 
specifically, the subscriber station obtains the measurements from one or more transmissions 
from base stations, sectors, GPS satellites, or any combination thereof, and provides them to the 
PDE, which then estimates the position of the subscriber station. 

[0023] Other systems, methods, features and advantages of the invention will be or will 
become apparent to one with skill in the art upon examination of the following figures and 
detailed description. It is intended that all such additional systems, methods, features and 
advantages be included within this description, be within the scope of the invention, and be 
protected by the accompanying claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0024] The components in the figures are not necessarily to scale, emphasis instead 
being placed upon illustrating the principles of the invention. In the figures, like reference 
numerals designate corresponding parts throughout the different views. 
[0025] Figure 1 is a simplified diagram of an example position location system. 
[0026] Figure 2 illustrates the process of determining a plurality of potential position 
estimates each from different measurement subsets. 

[0027] Figure 3 is a simplified diagram of an example position location system 
representing an example application of the subject invention. 

[0028] Figure 4 is a time domain representation of an example of a correlation function 
for a CDMA pilot signal. 

[0029] Figure 5 is a time domain representation of an example of a correlation function 
for a CDMA pilot signal which reflects the effects of a multi-path condition 
[0030] Figure 6 illustrates an example scenario which might give rise to the correlation 
function of Figure 5. 

[0031] Figure 7 is a time domain representation of an example of a correlation function 
for a CDMA pilot signal subject to a non line of sight condition. 

[0032] Figure 8 illustrates an example scenario which might give rise to the correlation 
function of Figure 7. 
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[0033] Figure 9 is a tints domain representation of an example of a correlation function 
for a CDMA pilot signal subject to multi-path and non line of sight conditions. 
[0034] Figure 10 illustrates an example scenario which might give rise to the correlation 
function of Figure 9. 

[0035] Figure 1 1 A is a simplified flowchart illustrating a first embodiment of a method 
according to the invention in which a position estimate is selected based on local topography. 
[0036] Figure 1 IB is a simplified flowchart illustrating a second embodiment of a method 
according to the invention in which a position estimate is selected responsive to refined costs 
determined based on local topography for a plurality of position estimates. 
[0037] Figure 12 is a simplified flowchart illustrating a third embodiment of a method 
according to the invention in which any unreliable measurements are excluded from 
consideration based on consideration of local topography. 

[0038] Figure 13 is a simplified flowchart illustrating a fourth embodiment of a method 
according to the invention in which a position estimate is selected responsive to reliability 
metrics determined based on consideration of local topography. 

[0039] Figure 14 is a flowchart illustrating an implementation example of a method 
according to the invention. 

[0040] Figure IS is a simplified block diagram of a subscriber station. 

DETAILED DESCRIPTION 

Example Application 

[0041] Referring to Figure 3, an example application of the invention is illustrated. In this 
example application, subscriber station 100 is a component of a wireless communication system 
such as but not limited to cellular, fixed wireless, PCS, and satellite communications systems. 
Moreover, the wireless communications system may provide for multiple access in accordance 
with CDMA, TDMA, FDMA, or GSM multiple access protocols, or combinations thereof. 
[0042] One or more base station(s) or sector(s) 102a, 102b, and 102c are employed in the 
wireless communications system. Each base station or sector 102a, 102b, 102c transmits a pilot 
signal which is modulated with a repeating pseudo-random noise (PN) code which uniquely 
identifies that base station or sector. For IS-95 compliant CDMA systems, the PN code is a 
sequence of 32,768 chips, which is repeated every 26.67 msec. 
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[0043] One or more GPS satellites 106a, 106b may also be visible to the subscriber 
station 100 or position determination entity (PDE) 104. Each of the GPS satellites also 
transmits a signal which is modulated with a repeating PN code which uniquely identifies that 
satellite. In current GPS systems, the PN code is a sequence of 1,023 chips, which is repeated 
every millisecond. 

[0044] A parameter estimator within subscriber station 100 is configured to estimate 
various parameters of the pilot signals transmitted from the base stations or sectors 102a, 102b, 
and 102c and/or the signals transmitted from the GPS satellites 106a, 106b. Such parameters 
may include TOA, time of transmission, energy per chip divided by total received power (EJIc)> 
and reliability metrics for any of the foregoing, including but not limited to root mean squared 
error (RMSE) for one or more of these parameters, such as the TOA estimate. 
[0045] The parameters, once estimated, and the reliability metrics if computed, are 
provided to PDE 104, which estimates the location of subscriber station 100 responsive thereto. 
(The PDE 104 may be a server in a public computer network such as the Internet or other 
TCP/IP network, a private network, or some other network.) As part of this process, the PDE 
104 may weight the parameters by their corresponding reliability metrics and estimate the 
location of subscriber station 100 using the weighted parameters. Once estimated, the position 
of the subscriber station 100 may be downloaded to it, or the position information may be 
transmitted to an emergency operator in the event of a 91 1 or other emergency call. 
[0046] The PDE 104 may estimate the location of the subscriber station 100 from 
measurements derived from transmissions from base stations or sectors 102a, 102b, and 102c, 
or, to increase accuracy, from the combined measurements of one or more of the base station(s) 
or sectors) 102a, 102b, 102c and one or more of the GPS satellites) 106a, 106b. 
[0047] The PDE 104 may provide other forms of assistance to the subscriber station 100. 
For example, PDE 104 may continuously track the GPS satellites, and provide assistance to the 
subscriber station 100 in locating the signals transmitted from them. This avoids the need to 
have the subscriber station 100 undergo time-consuming "cold start" procedures for locating the 
satellites when it is powered up. 

[0048] In one embodiment, the parameter estimator within subscriber station 100 
comprises correlation logic and analysis logic. For purposes of this disclosure, the term "logic" 
means hardware, software, or a combination of hardware and software. A composite signal 
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comprising a combination of multiple pilots signals transmitted from multiple base stations, 
sectors, and/or GPS satellites is received by the correlation logic. Each pilot signal is modulated 
with an identification code which, in one example, is a PN code. The identification code may 
be modulated onto the signal on a one-time or repeating basis. The identification code used to 
modulate a pilot signal uniquely identifies the source of the pilot signal, whether it be a base 
station, sector, or a GPS satellite. 

[0049] The correlation logic is configured to determine, using an integration time I, the 
correlation between the composite signal and shifted versions of an identification code. 
Obviously, the greatest degree of correlation will be exhibited if the identification code used 
by the correlation logic matches that modulated onto a pilot signal which is contained in the 
composite signal. The correlation logic outputs a correlation function which represents the 
correlation within a search window between the signal and shifted versions of the 
identification code. 

[0050] In one example, each sample of the composite signal S is a complex number 
having in-phase (I) and quadrature (Q) components, and the signal S comprises a plurality of 
pilot signals each modulated with a PN code. In one implementation, the correlation logic 
determines a correlation value C, which depends on the PN code and the shift s in the PN code 
which are used. This correlation value is coherent, i.e., retains phase information, and is a 
complex number which can be expressed as follows: 

C(PN,s)= £>tf(i-.y)«S(i) (4) 

i=*+0 

where N is the (coherent) integration time in terms of chips, S(i) are samples of the received 
signal, and k is an arbitrary start time. In this implementation, the integration time I is the 
coherent integration time N. 

[0051] In a second implementation, the correlation logic determines a correlation value 
C, which is a real number and is derived by non-coherently, i.e., not retaining phase 
information, combining M successive coherent integrations, each conducted over N chips. In 
this implementation, the correlation value C may be expressed as follows: 

C(JW,*) = £ £|PiV(«-*).5(/)| (5) 

In this implementation, the integration time I is the product of N and M. 
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[0052] The range of the shift s that is tested is the search window W. The W values 
C(PN, s) which result together form a correlation function F(PN, s), which represents the 
degree of correlation between the signal S and the shift s of the PN code (where the shift s is 
expressed in terms of chips) over a desired search window W. In the case where the PN code 
is repeatedly modulated onto the signal, the correlation function F(PN, s) will be periodic. 
[0053] Figure 4 illustrates an example of a correlation function F(PN, s) for a pilot 
signal in a CDMA wireless communication system. The window size (in terms of chips) in 
this example is 8, and the window is assumed to be centered at the origin 406. The 
horizontal axis 402 represents the shift of the PN code (expressed in terms of chips), and the 
vertical axis 404 represents the correlation function F(PN, s) (expressed in terms of 
Energy(dB)). As illustrated, the peak 408 of the function in this example is located at the 
origin 406. This correlation function may constitute a component pulse in a real world 
correlation function. 

[0054] If the peak 408 of the correlation function is sufficiently high and sufficiently 
distinguishable from noise and other distortions (to be discussed) that the risk of a false 
positive is reduced to an acceptable level, then a successful measurement has been obtained. 
The correlation logic may repeat this process with other PN codes until at least four 
measurements have been successfully obtained. 

[0055] For each measurement that has been successfully obtained, the analysis logic 
may then derive a propagation time estimate which, in one example, is the time/distance 
associated with the earliest non-sidelobe peak of the correlation function within the search 
window. The time/distance associated with the earliest non-sidelobe peak is generally used 
since it often corresponds to a Kne-of-sight transmission, and the peak associated with the 
line-of-sight transmission is usually needed for an accurate TOA estimate. The analysis logic 
may also derive a reliability metric such as RMS error for the TOA estimate, and weight the 
TOA estimate with a weight derived from the reliability metric. 

[0056] The analysis logic may provide PDE 104 with four or more weighted TOA 
estimates. Responsive to this information, PDE 104 estimates the location of subscriber 
station 100, and then downloads this information to the subscriber station 100. 
[0057] The subject invention may be employed within PDE 104 and utilized for the 
purpose of estimating the location of the subscriber station 100 responsive to the 
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measurements (weighted or unweighted) which have been provided to it. However, it should 
be appreciated that other applications are possible where the position determination function 
is performed directly within the subscriber station 100 or within another entity, such as one 
or more of the base stations or sectors 102a, 102b, 102c. In that event, the subject invention 
may be employed within any of these other entities within the wireless communications 
system. 

Background Regarding The Effects of Multi-path And/Or Lack of Line of Sight 
[0058] Figure S illustrates an example of a correlation function F(PN, s) for a pilot 
signal in a CDMA communications system where the effects of multi-path are depicted. The 
correlation function for a weaker earlier-in-time independent arrival, represented in the figure 
with a dashed line, is superimposed on the correlation function for a stronger independent 
arrival, which is represented in the figure with a solid line. The additive combination of these 
two correlation functions might represent the correlation function which is derived from a 
real world composite signal subject to the effects of multi-path or which is a multi-path 
signal. For purposes of this disclosure, an "independent" arrival includes one of multiple 
signals from different sources which arrive at a destination as well as one of multiple 
renderings of the same signal which arrive at a destination through different routes. 
[0059] Figure 6 illustrates a situation which might give rise to the correlation function 
depicted in Figure 5. As illustrated, the direct line-of-sight path between base station or 
sector 102 and subscriber station 100 is blocked by an obstruction 614 (a tree in the 
example), which attenuates but still allows passage of the pilot signal. At the same time, due 
to multi-path, the same pilot signal may reflect off of another obstruction 616 (a building in 
the example) and be received by the subscriber station 100 with less attenuation than the 
direct line-of-sight arrival. This reflection from obstruction 616 might give rise to the 
correlation function represented in Figure 5 with a solid line, while the direct line or sight 
rendering through obstruction 614 might give rise to the correlation function represented in 
Figure 5 with a dashed line. 

[0060] Turning back to Figure 5, the multi-path condition compounds the difficulty of 
locating the earliest non-sidelobe peak of the correlation function, which, in Figure 5, is the 
peak identified with numeral 512. The reason is that this peak may be difficult to distinguish 
from the sidelobes of the correlation function representing the reflected signal in Figure 6, 
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which sidelobes are identified with numeral 518. Accordingly, there will be a risk that the 
analysis logic in the parameter estimator will either erroneously interpret such a sidelobe as 
the peak of the correlation function for the line of sight arrival, or would erroneously 
interpret the peak of this correlation function as a sidelobe from a later in time arrival. Either 
way, the ensuing time-of-arrival estimate will be skewed. 

[0061] Figure 7 depicts an example of a correlation function derived from a pi lot signal 
in a CDMA wireless communications system which reflects the effects of a non-line of sight 
condition, such as might result from the scenario depicted in Figure 8. There, as illustrated, 
obstruction 804 completely blocks a line of sight transmission from base station 102 to 
subscriber station 100. However, a non-line of sight transmission from base station 102 
reflects off building 802 and is received at subscriber station 100. This non-line of sight 
transmission gives rise to the correlation function depicted in Figure 7. 
[0062] Here, an accurate estimate of time-of-arrival is difficult to obtain since a peak 
corresponding to a line of sight transmission is not present in the correlation function. 
Moreover, since peak 702 is the earliest non-sidelobe peak depicted in Figure 7, there is a 
risk that the analysis logic will use the time/distance associated with peak 702 in determining 
a time-of-arrival estimate for the pilot. Since this peak is not associated with a line of sight 
transmission, the time-of-arrival estimate will be skewed. 

[0063] Figure 9 illustrates a correlation function reflecting a combination of multi-path 
and line-of-sight conditions which might result in a CDMA wireless communications system. 
Figure 10 illustrates an example of a scenario which might give rise to the correlation 
function illustrated in Figure 9. As illustrated, a line of sight transmission from base station 
102 to subscriber station 100 is completely blocked by building 604. However, two reflected 
transmissions from base station 102 are received at subscriber station 100. The first, which is 
reflected off of building 612, is the earlier-in-time arrival and is assumed to give rise to the 
peak 902 illustrated in phantom in Figure 9. The second, which is reflected off of building 
602, is the later-in-time arrival and is assumed to give rise to the peak 904 illustrated in 
Figure 9. 

[0064] Here, an accurate estimate of time-of-arrival is again difficult to obtain since a 
peak corresponding to a line-of-sight transmission is not present in the correlation function. 
Moreover, since peak 902 is the earliest non-sidelobe peak, there is a risk the analysis logic 
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will use the time/distance associated with it in estimating time-of-arrival. Since this peak 
does not correspond to a line-of-sight transmission, the time of analysis estimate will 
accordingly be skewed. 
Embodiments of the Invention 

[0065] The invention provides a method of deriving a position estimate of a receiver from 
a plurality of parameter measurements. Each of the parameter measurements is derived from a 
correlation function, which in turn is derived by correlating a signal received at the receiver with 
an identification code uniquely identifying a transmitter visible to the receiver. 
[0066] Referring to Figure 11 A, a first embodiment of the method is illustrated, hi step 
1 100, this embodiment of the method comprises deriving a plurality of position estimates, each 
from a different subset of the plurality of measurements. From step 1 100, the method proceeds 
to step 1 102. There, one of the position estimates is selected based on local topography. 
[0067] A second embodiment of the method is illustrated in Figure 11B. As illustrated, 
this embodiment of the method begins with step 1104, which comprises deriving a plurality of 
position estimates, each from a different (although possibly overlapping) subset of the plurality 
of measurements which have been obtained. For example, if six measurements A, B, C, D, E, 
and F have been obtained, step 1104 might comprise deriving six position estimates from each 
the following measurement subsets: 1) A, B, C, D, E; 2) A, B, C, D, F; 3) A, B, C, E, F; 4) A, 
B, D, E, F; 5) A, C, D, E, F; 6) B, C, D, E, F. 

[0068] From step 1104, the method proceeds to step 1106, which comprises deriving the 
cost of each of the position estimates determined in step 1104. For purposes of this disclosure, 
the "cost" of an estimate is any measure of the consistency of the measurements used in 
deriving that estimate, and includes, but is not limited to, the RMS of the residual measurement 
errors. 

[0069] From step 1 106, the method proceeds to step 1 108, which comprises refining each 
of the costs determined in step 1 104 based on local topography. 

[0070] In one implementation, each cost is refined using the following process. First, one 
of the measurements used to derive the corresponding position estimate is selected, and a 
database representing the local topography is accessed to identify any obstacles in the path of 
the line of sight transmission from the corresponding reference source and the receiver. 
Optionally, any obstacles which may give rise to reflected transmissions at the receiver are also 
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identified. From this information, the expected correlation function for the reference source is 
derived and compared to the actual correlation function which is obtained. Based on this 
comparison, the cost for the estimate as determined in step 1106 might be adjusted. This 
process is repeated for each of the measurements used to derive the position estimate. The 
result is the refined cost for the position estimate. 

[0071] In one implementation example, the obstacles to a line of sight transmission are 
identified, and the expected amount of attenuation at the receiver from a line of sight 
transmission from the reference source is determined. Rom this information, the expected peak 
amplitude for the correlation function resulting from the line of sight transmission is then 
determined and compared with the amplitude of the earliest non-sidelobe peak of the actual 
correlation function. If there is a divergence, the cost may be adjusted upwards, with the 
amount of the adjustment depending on the degree of divergence. If there is a match or a near 
match, the cost may be adjusted downwards or kept the same. This process is repeated for each 
of the measurements from which the position estimate is derived. The result is the refined cost 
for the position estimate. 

[0072] Consider, for example, the situation where the actual peak amplitude exceeds the 
expected peak amplitude. Here, a multi-path condition could be present. Hence, the cost might 
be increased to reflect that. As another example, consider a situation where the actual peak 
amplitude is less than expected. Here, unexpected obstacles could be present in the line of sight 
path from the transmitter. Hence, the cost might be increased to reflect that. 
[0073] In a second implementation example, the expected shape of the correlation 
function from a line of sight transmission is estimated based on consideration of the obstacles 
expected to be present in a line of sight transmission. The estimated shape could also reflect the 
expected effects of multi-path based on consideration of obstacles expected to provide reflected 
transmissions to the receiver. Once determined, the estimated shape is compared with the actual 
shape. Again, the cost could be adjusted in response to the comparison. As before, this process 
would be repeated for each of the measurements from which the position estimate is derived. 
The result is the refined cost for the position estimate. 

[0074] Turning back to Figure 1 IB, from step 1 108, the method proceeds to step 1110. In 
step 1 110, a position estimate is selected responsive to refined costs. In one example, this step 
comprises selecting the position estimate with the lowest refined cost. 
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[0075] Figure 12 illustrates a third embodiment of the method of the invention. In this 
embodiment, the method begins with step 1204, which comprises excluding any unreliable 
measurements (from the plurality of measurements which have been obtained) based on 
consideration of local topography. 

[0076] In one implementation, this step comprises accessing, for each measurement, a 
database representing local topography to identify any obstacles expected to be present is a line 
of sight transmission from the corresponding transmitter and/or expected to give rise to a 
reflected transmission at the receiver. Responsive to this information, the expected correlation 
function for the transmitter is computed and compared with the actual correlation function. 
Responsive to this comparison, the measurement is either excluded as unreliable, or retained. 
[0077] In one implementation example, the expected attenuation of a line of sight 
transmission from the transmitter is computed, and the expected amplitude of the correlation 
function resulting from the line of sight transmission derived therefrom. This estimate is 
compared with the amplitude of the earliest non-si delobe peak, in the actual correlation function. 
If there is a significant enough divergence, the measurement may be rejected as unreliable. 
[0078] In a second implementation example, the expected shape of the correlation 
function for the transmitter is computed and compared to the shape of the actual correlation 
function. Based on this comparison, the measurement may be rejected as unreliable. 
[0079] Turning back to Figure 12, from step 1204, the method proceeds to step 1206. In 
step 1206, the method comprises deriving position estimates from different subsets of the 
remaining measurements. From step 1206, the method proceeds to step 1208, which comprises 
deriving costs for each of these position estimates. From step 1208, the method proceeds to step 
1210, which comprises selecting one of the estimates based on the costs. In one 
implementation, this step comprises selecting as the best estimate the position estimate with the 
lowest cost 

[0080] Figure 13 illustrates a fourth embodiment of the method of the invention. In this 
embodiment, the method begins with step 1304, which comprises deriving a plurality of 
position estimates, each from a different subset of the measurements which have been obtained. 
[0081] From step 1304, the method proceeds to step 1306, which comprises using local 
topography to compute reliability metrics for each of the position estimates derived in step 
1304. 
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[0082] From step 1306, the method proceeds to step 1308, which comprises selecting a 
position estimate as the best estimate responsive to the reliability metrics determined in step 
1304. 

[0083] The reliability metrics may represent a proportional or inversely proportional 
measure of the reliability of the corresponding position estimate. In the case where the 
reliability metric is directly proportional to the reliability of the position estimate, step 1308 may 
comprise selecting as the best estimate the position estimate with the largest reliability metric. 
In the case where the reliability metric is an inverse measure of the reliability of the position 
estimate, step 1308 may comprise selecting as the best estimate the position estimate with the 
smallest reliability metric. 

[0084] In one implementation, step 1306 comprises accessing, for each measurement used 
to derive a position estimate, a database representative of local topography to identify any 
obstacles in a line of sight transmission from the corresponding transmitter and/or any obstacles 
which could result in reflected transmissions at the receiver. Responsive to this information, the 
expected correlation function for the transmitter is obtained and compared with the actual 
correlation. Responsive to this comparison, the reliability metric for the estimate may be 
adjusted. This process is repeated for each of the measurements from which a position estimate 
is derived. The result is a reliability metric for the position estimate. In like manner, the 
reliability metric for each of other position estimates is derived. 

[0085] In one example, for each measurement used to compute a position estimate, the 
database is accessed to identify any obstacles in a line of sight transmission from the transmitter. 
From this information, the expected attenuation in a line of sight transmission from the 
corresponding transmitter is estimated, and the expected peak amplitude of the correlation 
function for the line of sight transmission determined. The expected peak amplitude is 
compared with the amplitude of the earliest non-sidelobe peak in the actual correlation function, 
and the reliability metric for the estimate adjusted accordingly. For example, if the comparison 
indicates a multi-path or non-line of sight condition is present or likely, the reliability metric 
could be adjusted upwards or downwards, depending upon if the metric is an inverse or 
proportional measure of reliability. In another example, for each measurement, the comparison 
involves comparing the expected and actual shape of the corresponding correlation function and 
adjusting the correlation function accordingly. 
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[0086] Again, this process is repeated for each of the measurements from which the 
position estimate is derived, thus yielding the reliability metric for the estimate. This overall 
process is then repeated for each of the estimates obtained in step 1304 to yield reliability 
indices for all the estimates. 

[0087] Figure 14 illustrates an implementation example of the method of the invention. 
[0088] The process starts with task 1400 and proceeds to task 1401 . 
[0089] In task 1401, a set of measurements is obtained. The elements of the measurement 
set include measurements from a number of reference points, for example GPS satellites, 
cellular base stations, cellular terrestrial measurements, etc. Elements of the measurement set 
may also represent multiple measurements for any given reference point This may involve 
multiple independent or correlated measurements, or multiple hypotheses for a single 
measurement, such as multiple correlation peaks. Each measurement can be provided with a 
weighting factor, such as an RMS error estimate. 

[0090] Rom task 1401 , the method proceeds to task 1402. Task 1402 queries whether all 
possible suitable subsets of the measurement set have been taken. A suitable subset would 
contain sufficient number of independent measurements to derive position solution(s). The 
solution does not need to be fully determined. A suitable subset typically does not contain 
measurements that are known to be inconsistent, such as multiple correlation peaks. 
[0091] If the answer is no, control passes to task 1403. In task 1403, a subset is selected. 
[0092] Task 1404 follows task 1403. In task 1404, a position solution is computed by 
using a suitable navigation algorithm. The navigation algorithm may use the measurement 
weighting factors in the computation. The output of the navigation algorithm may be a single 
point in the case of a fully determined or over-determined solution, or a finite set of points if 
ambiguous solutions exist, or even an infinite set of points in the case of an under-determined 
solution. An example for the latter is when pseudo-range measurements for three GPS satellites 
are used, in which case the solution will be a set of points forming a curve. Such sets of points 
should be delimited by using appropriate hypotheses. For example, a non-horizontal curve can 
be delimited by a certain height interval assumption (altitude aiding). The extent of the limit is 
determined by a trade-off between computational time and required reliability. 
[0093] Task 1405 follows task 1404. In task 1405, the set of points determined in task 
1404 is extended, as dictated by inherent measurement uncertainties. Each solution point is 
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typically extended to an uncertainty ellipsoid. The parameters describing the ellipsoid are 
results of the navigation algorithm, and are dependent on the individual measurement RMS 
estimates and the geometry of the reference points. If the navigation solution yielded more than 
one point, then the combined extension is determined as the union of the uncertainty areas 
assigned to every point. 

[0094] Task 1406 follows task 1405. Li task 1406, the extended set derived in task 1405 
is divided into regions. Each region is represented by its geometric center in the algorithm 
hereinafter. The size of the regions is governed by a trade-off between computational time and 
required accuracy. 

[0095] From task 1406, control passes to task 1407. In task 1407, a query is made 
whether all regions woe tested. 

[0096] If the answer is no, task 1408 is performed. In task 1408, a region is selected. 
[0097] Task 1409 follows task 1408. Task 1409 determines an a priori probability 
associated with the selected region. The a priori probability is computed based on: 

1. Information about the particular location pertaining to the region. This 
encompasses all a priori knowledge of the local environment. For example, if a region has had 
more frequent location requests associated with it in the past, it should have a higher assigned 
probability. Also, if an under-determined solution was used in task 1404, then, when delimiting 
the resulting infinite set of points, the resulting set typically doesn't have a uniform probability 
distribution. For example, if altitude aiding was used, then typically the probability is less at the 
extrema of the altitude range than at the middle. 

2. Location of the region relative to the navigation solution. When extending the 
solution in task 1405, the uncertainty region has a certain probability distribution associated 
with it (for example, a 3-D Gaussian distribution for an uncertainty ellipsoid), which can be 
used as a factor in the a priori probability of the region. If the navigation solution yielded more 
than one point, then a particular region can be in the uncertainty region of a number of solution 
points. In this case, the combined probability is computed based on conditional probability 
mles. 

3. The sum of the weights of the measurements used in the navigation algorithm 

4. The cost given by the navigation algorithm for the solution 
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[0098] Task 1410 follows task 1409. In task 1410, a query is made whether all 
measurements of the measurement subset chosen in task 1403 have been considered. 
[0099] If no, task 141 1 is performed. In task 141 1, a measurement is selected. 
[0100] Task 1412 follows task 1411. In task 1412, a predicted receive signal strength is 
calculated. The calculation is based on the knowledge of the receiver location (which is the 
center of the region selected in task 1408), the location of the signal source (can be determined, 
for example, from the Ephemeris and possibly differential corrections for GPS satellites), and 
also on the knowledge of local topography. After drawing a line of sight from the receiver 
location to the signal source location, any possible obstacles can be identified. Based on the 
type of the obstacle and the length of the path blockage, an obliquity factor is determined and 
incorporated in an expected attenuation value. Knowledge of local topography can include 
knowledge of architectural objects, as well. For example, if the path is blocked by a building, 
great signal attenuation should be expected. The attenuation can be modeled as highly 'non- 
linear'. For example, if the required penetration depth of the signal inside the building is small, 
then a receiver close to a window, receiving the signal through the window can be assumed, 
which would mean small signal loss. The same receiver would be expected to see a much 
greater signal loss if the signal source's location is diagonally opposite, because then the signal 
would have to propagate through internal building structure. These assumptions can be easily 
captured in a 'non-linear' path-loss model, where a long path blockage would have an even 
greater expected attenuation compared to a short path blockage than what a homogeneous 
model would suggest. Along with the expected attenuation, an estimation error RMS estimate is 
calculated. If the path crosses multiple objects, then a combined signal attenuation and 
combined RMS is calculated. 

[0101] The attenuation also includes free space propagation loss, which can be calculated 
since the assumed recei ver-to-signa] source distance is known. 

[0102] Note that since the receiver typically sees cellular base stations at low elevation 
angles, it may be beneficial to approximate the attenuation for cellular base stations by simply 
using a generic attenuation vs. receiver-to-signal source distance function. The function can be 
chosen based on knowledge of the environment or based on collected statistical data. For 
satellites that are seen typically at higher elevations, the particular path and the expected 
blockage should be used in the calculations. 



WO 03/089954 



PCT/DS03/11475 



20 

[0103] Once the expected attenuation is determined, it is multiplied by the nominal source 
power to obtain the expected receive signal strength. For satellites the nominal source power is 
known, for cellular base stations, it can possibly be read from overhead messages broadcast by 
the base station. In the absence of such a broadcast, an estimate can be used. Any uncertainty 
in the source power should be reflected in the receive signal strength RMS error estimate. 
[0104] From task 1412, control passes to task 1413. In task 1413, the actual receive 
signal strength and the predicted receive signal strength is compared, and a conditional 
probability is computed. A low probability could be the result of either seeing a signal with a 
much greater power than expected, which would indicate that the detected signal is probably a 
multi-path (since the line-of-sight signal is expected to be blocked); or seeing a signal with 
much less power than expected. In either case, the confidence in the region selected in task 
1408 (and being the location where the measurement selected in task 141 1 is taken) is low. 
[0105] If the query in task 1410 indicates that all measurements in the subset selected in 
task 1403 have been considered, a combined probability is computed in task 1414. This 
computation is based on the results obtained in task 1413. 

[0106] From task 1414, control passes to task 1415. In task 1415, a combined probability 
is computed, based on the results of tasks 1409 and 1414, and control then passes back to task 
1407. 

[0107] Once the query in task 1407 indicates that all regions have been tested, task 1416 
is performed. There, the region for which the maximum probability was computed in task 1415 
is selected and recorded. The geometric center of the recorded region will be the location 
solution associated with the measurement subset selected in task 1403, and the probability of the 
region will be the probability assigned to the measurement subset. From task 1416, control 
passes back to task 1402. 

[0108] If the query in step 1402 indicates that all measurement subsets have been tested, 
task 1417 is performed. There, the measurement subset with the highest probability is selected, 
and the associated location is declared the position solution. 
[0109] The process then exits at task 1418. 

[0110] Not all steps described above are required for operation. If certain sets of 
information are not available, then the corresponding steps may be omitted. 
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[0111] Any of the foregoing methods may be tangibly embodied in a variety of forms, 
including but not limited to, a form where a series of instructions embodying the method is 
stored on or in a memory or a server in a computer network such as the Internet, where the 
method is embodied as synthesized logic, or where the method is embodied as a computer 
program product, i.e., a code segment or module. 

[0112] Furthermore, systems are possible comprising a processor configured to access a 
memory embodying any of the foregoing methods to obtain and execute the instructions stored 
thereon. Such a system may be used in a variety of environments or applications. For example, 
in one application, such a system is situated in a subscriber station and is configured to locate 
the position of the subscriber station from one or more transmissions from base stations, sectors, 
GSP satellites, or any combination thereof. In another application, such a system is situated in a 
position determination entity (PDE), and is configured to locate the position of a subscriber 
station which obtains measurements from one or more transmissions from base stations, sectors, 
GPS satellites, or any combination thereof, and provides them to the PDE, which then 
determines an estimate of the position of the subscriber station. 

[0113] For purposes of this disclosure, a "processor" may be any device capable of 
executing a series of instructions embodying a process, including but not limited to a computer, 
microprocessor, an ASIC, finite state machine, DSP, or some other mechanism. 
[0114] Moreover, as used herein, a "memory" may be any device readable by a 
processor and capable of storing a series of instructions embodying a process, including but 
not limited to RAM, ROM, EPROM, EEPROM, PROM, disk (hard or floppy), CD-ROM, 
DVD, flash memory, etc. 

[0115] Referring to Figure 15, a simplified block diagram of one example implementation 
of a subscriber station is illustrated, hi this particular example, numeral 1502 identifies a 
processor configured to execute software instructions, and numeral 1504 identifies a memory 
for holding software instructions and data, including software instructions embodying the 
method of the invention. 

[0116] Numeral 1506 is persistent storage for holding provisioning information useful for 
roaming and/or system registration and de-registration. 

[0117] Numerals 1508 and 1510 identify, respectively, a keypad and display, and are 
typically provided as part of a user interface. 
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[0118] Numerals 1512 and 1514 identify, respectively, a microphone and speaker, and are 
typically provided to support uses involving transceiving voice. 

[0119] Numeral 1516 identifies a radio capable of receiving and transmitting information 
at RF frequencies over a wireless interface according to one or more multiple access protocols, 
such as but not limited to CDMA, GSM, WCDMA, AMPS, TDMA, PCS, or combinations 
thereof. 

[0120] Numeral 1518 identifies a modem, which is capable of upconverting voice (and 
data) at baseband frequencies to RF frequencies, and downc on verting RF frequencies to 
baseband frequencies. 

[0121] Numeral 1520 identifies a correlator, which is configured to derive a correlation 
function from a received signal in the manner previously described. 

[0122] Numeral 1522 identifies a database representative of local topography, which is 
accessible by processor 1502. 

[0123] In this particular example, measurements are derived by processor 1502 from the 
correlation functions provided by correlator 1520. Processor 1502 then executes the software 
instructions stored in memory 1504, while accessing database 1522, to perform the method of 
the invention. At the conclusion of this process, an estimate of the position of the subscriber 
station is determined. 

[0124] When the technique represented by any of these embodiments, implementations or 
examples is applied to measurements or position estimates derived from such measurements, 
performance is greatly improved since distortions caused by factors such as multi-path and lack 
of line of sight may now be accurately accounted for. 

[0125] While various embodiments of the invention have been described, it will be 
apparent to those of ordinary skill in the art that many more embodiments and 
implementations are possible that are within the scope of this invention. 



What is claimed is: 
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CLAIMS 

1 . A method of deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

deriving a plurality of position estimates, each from a different subset of the plurality 
of measurements; 

deriving a cost of each of the position estimates; 

refining the costs based on local topography; and 

selecting one of the position estimates responsive to the refined costs. 

2. The method of claim 1 wherein the parameter measurements are time-of- 
arrival measurements. 

3. The method of claim 1 wherein the identification codes are pseudo noise (PN) 

codes. 

4. The method of claim 1 wherein the cost of an estimate is the RMS of the 
residual errors of the measurements used in deriving that position estimate. 

5. The method of claim I wherein the refining step comprises performing the 
following substeps for each measurement: 

estimating, based on local topography, the expected correlation function for the 
corresponding transmitter; 

comparing the expected correlation function with the actual correlation function 
which is obtained at the receiver, and 

adjusting, based on the comparison, the cost for each of the position estimates derived 
from the measurement, if the divergence between the expected and actual correlation 
function is sufficiently great. 

6. The method of claim 5 wherein the comparing substep comprises comparing 
the expected peak amplitude with the amplitude of the earliest non-sidelobe peak of the 
actual correlation function. 
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7. The method of claim 5 wherein the comparing substep comprises comparing 
the expected shape of the correlation function with the actual shape of the correlation 
function. 

8. The method of claim 1 tangibly embodied as a series of instructions stored in 
a memory. 

9. A system comprising a processor and the memory of claim 8, wherein the 
memory is accessible by the processor, and the processor is configured to access and execute 
the instructions stored in the memory. 

10. A method of deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

excluding, based on local topography, any of the plurality of measurements which are 
determined to be unreliable; 

deriving a plurality of position estimates, each from a different subset of the 
remaining measurements; 

deriving a cost of each of the position estimates; and 

selecting one of the position estimates responsive to the costs. 

1 1 . The method of claim 10 wherein the identification codes are pseudo noise 
(PN) codes. 

12. The method of claim 10 wherein the measurements are time-of-arrival 
measurements. 

13. The method of claim 10 wherein the excluding step comprises performing, for 
each measurement, the following substeps: 

estimating, based on local topography, the expected correlation function for the 
corresponding transmitter; 

comparing the expected correlation function with the actual correlation function from 
which the measurement is derived; 

determining, based on the comparison, whether the measurement is unreliable; and 

excluding the measurement from further consideration if determined to be unreliable. 
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14. The method of claim 13 wherein the comparing substep comprises comparing 
the expected peak amplitude of the correlation function with the amplitude of the earliest 
non-sidelobe peak of the actual correlation function. 

15. The method of claim 13 wherein the comparing substep comprises comparing 
the expected shape of the correlation function with the actual shape of the correlation 
function from which the measurement is derived. 

16. The method of claim 10 tangibly embodied as a series of instructions stored in 
a memory. 

17. A system comprising a processor and the memory of claim 16, wherein the 
memory is accessible by the processor, and the processor is configured to access and execute 
the instructions stored in the memory. 

18. The method of claim 10 wherein at least one of the transmitters is a GPS 
satellite. 

19. The method of claim 10 wherein at least one of the transmitters is a base 
station in a wireless communications system. 

20. The method of claim 10 wherein at least one of the transmitters is a base 
station sector in a wireless communications system. 

21 . The method of claim 10 wherein the signal is received at the receiver through 
wireless transmission. 

22. The method of claim 21 wherein the wireless transmission occurs in 
accordance with a multiple access protocol. 

23. The method of claim 22 wherein the multiple access protocol is CDMA. 

24. The method of claim 22 wherein the multiple access protocol is GSM. 

25. The method of claim 22 wherein the multiple access protocol is WCDMA. 

26. A method of deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

deriving a plurality of position estimates, each from a different subset of the plurality 
of measurements; 
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deriving, based on local topography, a reliability metric for each of position 
estimates; and 

selecting, based on the reliability metrics, one of the position estimates. 

27. The method of claim 26 wherein the step of deriving a reliability metric for a 
position estimate comprises performing the following substeps: 

determining the expected correlation function corresponding to each of the 
measurements from which the position estimate is derived; 

comparing, for each of the measurements, the expected correlation function with the 
actual correlation function from which the measurement is derived; and 

deriving, based on the comparisons for each of the measurements, the reliability 
metric for the position estimate. 

28. The method of claim 27 wherein the comparing step comprises comparing 
expected and actual peak amplitudes. 

29. The method of claim 28 wherein the comparing step comprises comparing an 
expected peak amplitude for a line of sight transmission with the actual amplitude of the 
earliest non-sidelobe peak. 

30. The method of claim 27 wherein the comparing step comprises comparing 
expected and actual shapes of correlation functions. 

3 1 . The method of claim 26 tangibly embodies a series of instructions stored in a 
memory. 

32. A system comprising a processor and the memory of claim 31, wherein the 
memory is accessible by the processor, and the processor is configured to access and execute 
the instructions stored thereon. 

33. The system of claim 32 in a subscriber station. 

34. The system of claim 32 in a PDE. 

35. A method of deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

deriving a plurality of position estimates, each from a different subset of the plurality 
of measurements; and 
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selecting, based on local topography, one of the position estimates. 

36. A method of deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

a step for deriving a plurality of position estimates, each from a different subset of the 
plurality of measurements; and 

a step for selecting, based on local topography, one of the position estimates. 

37. A system for deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

first logic for deriving, from the plurality of measurements, a plurality of position 
estimates, each derived from a different subset of the measurements; 

second logic for deriving a cost of each of the position estimates; 

third logic for refining the costs based on local topography; and 

fourth logic for selecting one of the position estimates responsive to the refined costs. 

38. The system of claim 37 wherein the parameter measurements are time-of- 
arrival measurements. 

39. The system of claim 37 wherein the identification codes are pseudo noise 
(PN) codes. 

40. The system of claim 37 wherein the cost of an estimate is the RMS of the 
residua] errors of the measurements used in deriving that position estimate. 

41 . The system of claim 37 wherein the third logic is configured, for each 
measurement from which a position estimate is derived, to (a) estimate, based on local 
topography, the expected correlation function for the corresponding transmitter; (b) compare 
the estimated correlation function with the actual correlation function; and (c) adjust, based 
on the comparison, the cost for each of the position estimates derived from the measurement 
if the divergence between the expected and actual correlation functions is sufficiently great 
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42. The system of claim 41 wherein the third logic is configured to compare the 
expected peak amplitude of the correlation function with the amplitude of the earliest non- 
sidelobe peak of the actual correlation function. 

43. The system of claim 41 wherein the third logic is configured to compare the 
expected shape of the correlation function with the actual shape of the correlation function. 

44. A system for deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in rum is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

first logic for excluding, based on local topography, any of the plurality of 
measurements which are determined to be unreliable; 

second logic for deriving a plurality of position estimates, each from a different 
subset of the remaining measurements; 

third logic for deriving a cost of each of the position estimates; and 

fourth logic for selecting one of the position estimates responsive to the costs. 

45. The system of claim 44 wherein the identification codes are pseudo noise 
(PN) codes. 

46. The system of claim 44 wherein the measurements are time-of-arrival 
measurements. 

47. The system of claim 44 wherein the first logic is configured, for each 
measurement, to (a) estimate, based on local topography, the expected correlation function 
corresponding to the transmitter from which the measurement is derived; (b) compare the 
expected correlation function with the actual correlation function from which the 
measurement is derived; (c) determine, based on the comparison, whether the measurement is 
unreliable; and (d) exclude the measurement from further consideration if determined to be 
unreliable. 

48. The system of claim 47 wherein the first logic is configured to compare the 
expected peak amplitude with the amplitude of the earliest non-sidelobe peak of the actual 
correlation function from which the measurement is derived. 
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49. The system of claim 47 wherein the first logic is configured to compare the 
expected shape with the actual shape of the correlation function from which the measurement 
is derived. 

50. A system for deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

first logic for deriving a plurality of position estimates, each from a different subset of 
the measurements; 

second logic for deriving, based on local topography, a reliability metric for each of 
position estimates; and 

third logic for selecting, based on the reliability metrics, one of the position estimates. 

51. The system of claim 50 wherein the second logic is configured, for each 
position estimate, to (a) determine the expected correlation function corresponding to each of 
the measurements from which the position estimate is derived; (b) compare, for each of the 
measurements, the expected correlation function with the actual correlation function from 
which the measurement is derived; and (c) derive, based on the comparisons for each of the 
measurements, the reliability metric for the position estimate. 

52. The system of claim 51 wherein the second logic is configured to compare 
expected and actual peak amplitudes. 

53. The system of claim 52 wherein the second logic is configured to compare the 
expected peak amplitude for a line of sight transmission with the actual amplitude of the 
earliest non-sidelobe peak. 

54. The system of claim 51 wherein the second logic is configured to compare 
expected and actual shapes of correlation functions. 

55. The system of claim 50 in a subscriber station. 

56. The system of claim 50 in a PDE. 

57. A system for deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 
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first logic for deriving a plurality of position estimates, each from a different subset of 
the measurements; and 

second logic for selecting, based on local topography, one of the position estimates. 

58. A system for deriving a position estimate of a receiver from a plurality of 
parameter measurements, each derived from a correlation function which in turn is derived 
by correlating a signal received at the receiver with an identification code uniquely 
identifying a transmitter visible to the receiver, comprising: 

means for deriving a plurality of position estimates, each from a different subset of 
the plurality of measurements; and 

means for selecting, based on local topography, one of the position estimates. 
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